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The energies of the excited states of the Nucleon, A and Q. are computed in lattice QCD, using 
two light quarks and one strange quark on anisotropic lattices. The calculations are performed at 
three values of the pion mass: m K = 392(4), 438(3) and 521(3) MeV. We employ the variational 
method with a basis of about ten interpolating operators enabling six energies to be distinguished 
clearly in each irreducible representation of the octahedral group. We compare our calculations of 
nucleon excited states with the low-lying experimental spectrum. There is reasonable agreement 
for the pattern of states. 
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1. Introduction 

The goal of determining the spectrum of hadron masses from lattice QCD is addressed in 
this work by calculations of the excited state spectrum of the nucleon, A and Q. using anisotropic 
lattices.[l] Earlier excited-baryon analyses were based on quenched QCD [2] and two-light-flavor 
(Nf = 2) QCD [3]. In this work we use ensembles of gauge configurations developed in Ref. [4] for 
Nf = 2 + 1 QCD with two dynamical light quarks and one strange quark. The lattices are 16 3 x 128 
with spatial and temporal lattice spacings a s = 0.122 fm and a, = 0.035 fm. 

For families of particles with given isospin and strangeness, spectra are calculated in the six 
double-valued irreducible representations (irreps) of the octahedral group. There are three irreps 
for even-parity that are labeled with a g subscript (gerade) and three for odd-parity that are labeled 
with a u subscript (ungerade). They are: G\ g ,H g ,G2g-,G\ m H u and G2„. 

Sets of seven to eleven three-quark operators are used in each irrep and the variational method [5, 
6] is used to extract energies of six states. Most operators incorporate gauge-covariant displace- 
ments of the quarks relative to one another in order to obtain nontrivial shapes. [7] The recently 
developed "distillation" method [8] is used for quark smearing. We start with a large set of opera- 
tors in each irrep and then "prune" them to sets of about 10 that have the lowest condition numbers. 
That yields sets of approximately linearly-independent operators that are suitable for calculations 
based on diagonalizing a matrix of correlation functions. 

2. Results 

A detailed presentation of all of our results is given in Ref. [1]. Here we present selected 
results for the nucleon and A excited states. 

Plots of the nucleon effective energies, calculated as 

EM= l -,n(l^l), (2.1) 

where X(t) is an eigenvalue of the generalized eigenvalue problem, are shown in Figure 1 for the 
G\g and G\ u irreps. These plots show the values of E e s obtained from Eq. (2.1) as vertical bars and 
E e ff calculated using the fit function, 

X fit (t) = (1 -A)e- E{t -' o) +Ae- E ' {t -' \ (2.2) 

in place of A (?) in Eq. (2. 1) as curved dashed lines. Comparison of the curved dashed lines with the 
bars from the lattice ensembles shows the usefulness of two-exponential fits. The term Ae~ E 
models the contributions of higher energy states at early times allowing the exponential term ( 1 — 
A)e~ E ( t ~ t °* > to be determined over a larger fit window (*;,?/) than would be possible using a single 
exponential. Fit energy E and uncertainty of the fit energy, a, are shown by dashed horizontal 
lines at E + a and E — a extending over the fit window. Note that the statistics allow credible 
determinations of six energy levels in each irrep. This provides evidence that quark smearing 
based on "distillation" is effective with regard to suppressing high-frequency fluctuations in the 
gauge ensembles. 
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The energies obtained from the G\ g and G\ u effective mass plots of Fig. 1 are shown as boxes 
extending from E — a to E + a in Fig. 2. We show nucleon energies that are obtained in the 
same manner as shown in Fig. 1 for all irreps of the octahedral group and three pion masses. 
Experimental spectra are shown to the left of lattice energies for the spins and parities that have 
subductions to the lattice irreps. Lattice and experimental spectra are shown in Fig. 3 for the A 
family. See Ref. [1] for the £2 spectra. 

In the nucleon spectra, there is good evidence for a spin f state. We find nearly degenerate 
H u and Giu partner states, which is the signature of spin | . However, other spins are difficult 
to identify because there are many nearly-degenerate states, within uncertainties. It is a near-term 
goal within the collaboration to address spin identification by using operators that are subduced 
from continuum spins. 

Our lattice spectra show scant evidence for multiparticle states even though many energies lie 
above the relevant thresholds. This is probably because single-hadron operators are used. There 
is an inference for a multiparticle contribution in that we find four low-lying states in H u while 
there are three low-lying experimental states that have subductions to H u . There is a threshold for a 
multiparticle state in the same energy range. Our lattice results agree with the experimental pattern 
if one of the four low-lying H u states is multiparticle. However, we cannot identify the multiparticle 
states in the spectrum. It is a near-term goal within the collaboration to incorporate multiparticle 
operators that couple directly to such states. 

Some lattice states appear to be "squeezed" by the small lattice volume used. They show up at 
higher energies than would be the case in a larger volume. The G2 states require partner states in 
other irreps, such as H, in order to realize all the magnetic substates for a given spin. The partners 
should be close to the same energy. However, in the A spectra of Fig. 3 we find G2 states at high 
energies without suitable partners being evident. Possibly they have been "squeezed". It is also a 
goal to perform calculations of spectra at larger volumes. 

Although we do not attempt to extrapolate energies to m n = 140 MeV, it is evident from Figs. 2 
and 3 that the lowest-energy states on the lattice tend toward the energies of the physical resonances 
as the pion mass decreases. Decreasing the pion mass is an obvious goal but we recognize that it 
entails a more complex analysis for excited states that can decay. 

3. Summary 

This work represents a milestone in our long-term research program aimed at determining the 
spectra of baryons in QCD. It provides the first spectrum for N, A and Q. baryons based on Nf = 
2 + 1 QCD with high statistics. A large number of baryon operators is used to calculate matrices of 
correlation functions. They are analyzed using the variational method with fixed eigenvectors. The 
analysis provides spectra at three pion masses: m % = 392(4) MeV, 438(3) MeV and 521(3) MeV. 

The lattice volume and pion masses used give considerably higher energies than the experi- 
mental resonance energies. However, there is reasonable agreement of the overall pattern of lattice 
and experimental states. One exception is that almost all A states in the G2 irrep are too high. That 
may be caused by a volume that is too small for highly excited states. 
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Figure 1: Nucleon Gi g effective energies are shown for the lowest states in the upper six graphs. The 
effective energy increases from left to right along the first row and continues to increase from left to right 
along the second row. The lower six graphs show nucleon Gu, effective energies increasing in the same 
pattern. Calculations are for m K = 392(4) MeV. Vertical bars show the effective energy and the curved 
dashed line shows the effective energy calculated from the fit function. Horizontal dashed lines show the fit 
results for E ± a and their extent shows the fitting interval (ti,tf). 



5 



Nucleon, Delta and Omega excited state spectra 



Stephen J. Wallace 



> 




_ 2000 
> 



392 438 521 



> 




> 

<u 



392 438 521 




392 438 521 



> 




> 



392 438 521 




392 438 521 



Figure 2: Spectra for isospin | (nucleon family) at three values of m K in each irrep of the cubic group are 
compared with experimental spectra. Columns labeled by m n = 392, 438 and 521 MeV show lattice spectra. 
The Gi g and G\ u spectra in the m K = 392 MeV column are obtained from the plots of Fig. 1. Boxes extend 
from E — a to E + a. Two, three and four-star experimental resonances are shown to the left of lattice spectra 
in columns labeled by their J p values. Each J p value listed has a subduction to the lattice irrep shown. Each 
box for an experimental resonance has height equal to the full decay width and an inner box (color aqua) 
showing the uncertainty in the Breit-Wigner energy. Open boxes show the thresholds for multiparticle states. 
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Figure 3: Spectra for isospin I (A family) at three values of m K in each irrep of the cubic group are compared 
with experimental spectra. Columns labeled by m K = 392, 438 and 521 MeV show lattice spectra. Boxes 
extend from E — a to E + a. Two, three and four-star experimental resonances are shown to the left of 
lattice spectra in columns labeled by their J p values. Each J p value listed has a subduction to the lattice irrep 
shown. Each box for an experimental resonance has height equal to the full decay width and an inner box 
(color aqua) showing the uncertainty in the Breit-Wigner energy. 
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